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This is a review of applications of the Color Glass Condensate to the phenomenology of 
relativistic heavy ion collisions. The initial stages of the collision can be understood in 
terms of the nonperturbatively strong nonlinear glasma color fields. Wc discuss how the 
CGC framework can and has been used to compute properties of the initial conditions of 
AA collisions. In particular this has led to recent progress in understanding multiparticle 
correlations, which can provide a directly observable signal of the properties of the initial 
stage of the collision process. 
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1. Introduction 

The Relativistic Heavy Ion colhder (RHIC) at Brookhaven National Laboratory 
has been taking data from nucleus-nucleus, deuteron-nucleus and proton-proton 
collisions starting from 2000. This successful program has provided a wealth of in- 
formation on the properties of QCD matter at high energy density. For a review 
of the early experimental results we refer the reader to the experimental "white 
papers"^ by the RHIC collaborations. Since the earliest RHIC observations it has 
become clear that the produced matter is deconfined and cannot be understood in 
terms of phenomcnological models of low energy hadron physics. It is obvious that 
this will be even more true for heavy ion collisions at the LHC. On the other hand, 
there are many signs that the matter interacts far too strongly to be described 
by naive perturbative calculations. The first principles method of studying nonper- 
turbative QCD is to directly discretize the path integral on a lattice and evaluate 
it numerically. This is indeed the method of choice for the description of a static 
system, such as studying the equation of state. The lattice, however, is ill suited 
for understanding time dependent phenomena, and other methods must be devel- 
oped. One is then faced with a choice between two views, depending on whether 
"strongly interacting" in the context of gauge theory implies a large value of the 
coupling constant as or not. One option is that one should consider as to be large. 
Performing practical computations in this strong coupling limit has now become 
possible at least in some nonabelian gauge theories, if not yet in QCD, due to a 
duality transformation that makes it possible to translate the problem into solving 
a classical wave equation in 5 dimensional gravitjEI. There has been much activity 
in this direction recently, and gauge / gravity duality seems a promising tool to gain 
qualitative understanding of the properties of the medium in the plasma phase. 

A stong coupling theory fails, however, in describing some basic features of 
high energy hadronic collisions that characterize the initial stage of a heavy ion 
collision. The predominance of forward scattering in the angular distribution (i.e. 
the approximate boost invariance of particle production, the small baryon stopping 
in hadronic collisions, the existence of jets) point towards an interaction that is well 
described by a weak coupling gauge theory (see e.g. the discussion ii|2l). This leads 
to the point of view that we will assume in this paper; namely that for the initial 
conditions of heavy ion collisions the system can be described assuming that as 
is small. In spite of the coupling constant being small the system can be strongly 
interacting because it is dense; gluonic states have high occupation numbers and 
the interactions between them are highly nonlinear. One possible application of this 
idea is to assume the existence of a geometrical condition that restricts the number 
of particles in the central rapidity region when they begin to overlap in phase 
space. In a practical phenomcnological application one can then first calculate the 
initial pruduction of partons using a framework developed for a generic hadronic 
collision and then supplement this calculation with the geometric ("final state") 
saturation criterion. Examples of this line of reasoning are provided e.g. by the 
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EKRT niode3, which uses the usual cohinear pQCD factorization, or the parton 
percolation model^, based on string-like phenomenology. The line of thought that 
we will pursue here, however, is that the nonlinear interactions that control the 
behavior of the bulk of particle production are present already in the initial nuclear 
wavefunctiorl^l and lead to the emergence of an ("initial state") saturation scale; a 
transverse momentum scale that then manifests itself in particle production in the 
collision. 

In collisions of protons and nuclei the typical values of x that are probed in 
the wavefunction are x ~ px/^/s, where px is s typical transverse momentum scale 
of the particles being produced. Let us consider, in the center-of-mass frame of 
the collision, the nucleus moving in the +z direction. The parton with momentum 
fraction x will have longitudinal momentum p"*" ~ xy^/A and will thus probe the 
other, leftmoving, nucleus at a length scale Ax^ ^ A/{xy/s). The longitudinal size 
of the leftmover is Lorentz-contracted from A^^^Rp to ~ A^/'^Rp{Am]y /y/s). 

We see that if x <C l/{A^^^RpmN), the partons in the rightmoving nucleus will not 
be able to resolve the individual nucleons of the leftmoving one. The whole nucleus 
must therefore be treated as one coherent target, not as a collection of independent 
nucleons. 

The observation that the large x localized, valence-like, degrees of freedom are 
not resolved in the collision, but only the smaller x partons that they radiate, 
naturally leads to the idea of treating the two separately in an effective field theory 
approach. This effective field theory is known as the Color Glass Condensate (CGC) 
^ ^ (scc^ for reviews ancff^for a summary of the case for the CGC based on RHIC 
data.). The CGC describes a high energy hadron in terms of a color field (the small 
X gluons) radiated by an effective color current (the large x degrees of freedom). At 
high energy (small x) the radiation is enhanced by large logarithms of the energy 
and the number of gluons is large. When the occupation numbers of gluonic states 
in the wavefunction becomes large enough, ^ l/cts, the field can be described as a 
classical one, radiated by classical effective color charges. The classical color charges 
are stochastic random variables with a probability distribution Both the 

properties of the initial stage of a heavy ion collision and observables in DIS at small 
X can be computed in terms of these same classical gluon fields. The color charge 
distribution depends on nonperturbative input and cannot completely be computed 
from first principles. Its dependence on the energy scale (rapidity) that separates 
the large and small x degrees of freedom can, however, be computed and expressed 
in terms of a renormalization group equation. The distribution of color charges is 
a universal object; it can be measured in one process (ideally DIS) and then used 
as an independent input to make prediction for another one (say, the initial field 
configurations in a heavy ion collision). In this sense the situation is analogous 
to collinear factorized perturbation theory; there is a universal, nonperturbative 
distribution (color charge distribution or parton distribution function) , a separation 
scale (rapidity or virtuality) and a renormalization group equation derived from 
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first principles that describes the dependence of the nonperturbative input on this 
separation scale. 

2. Glasma initial state 

We now turn to the basic features of the glasmgllH fields in the initial stages of the 
collision. We shall show hot their structure as initially longitudinal boost invariant 
fields follows directly from the CGC description of the nuclear wavefunction. 

2.1. Description of a nucleus in the CGC 

In the CGC picture the small x degrees of freedom that contribute to bulk gluon 
production are treated as a classical color field radiated by (classical) color sources. 
As classical fields they obey the Yang-Mills equation of motion 

(1) 

with the source currents 

^6^+pi{^T,x-) + S'^-p2{^T,x+) (2) 

representing the large x part of the wavefunctions. The color charge densities 
are static in the sense that they do not depend on the relevant light cone time 
{x^ for pi^2 respectively). The color charges of the two nuclei are naturally in- 
dependent of each other. They are also localized on their respective light cones: 
Pi,2(xt,s^) ~ pi,2{'^t)S{x^). a naive way to understand this dependence is to 
argue that at high energy the nuclei are Lorentz-contracted to infinitesimal sheets 
in an interval Ax~ ^ Vv^- This argument does not, however, take into account 
the quantum nature of the color charged integrated into the effective description 
in terms of the p's. The color charge densities represent degrees of freedom that 
have a higher x, i.e. a higher longitudinal momentum than the ones described 
as the classical field. Consequently they are better localized in x~ (the conjugate 
variable to p"*") than the classical field and are, in the high energy (multi-Regge) 
kinematics we are working in, seen by the glasma field as infinitesimal sheets. The 
actual scale of this delta function can be estimated as follows. The longitudinal (and 
transverse) momenta of the gluons in the classical field ar^ ~ Qg. The evolution 
speed in rapidity is proportional to as; we can therefore assume that the sources 
are separated from the classical field by Ay ^ l/c^s units in rapidity. Thus they 
correspond to longitudinal momenta p~^ ~ e^^^'Qs, i.e. are localized in an interval 
Ax- ~ e-i/"=/(3s- This is the actual meaning of the delta function approximation 
used when writing pi{iiT, x^) ^ pi{xt)S{x^). As we shall see in the following, it is 
nevertheless important to initially maintain an explicit ^"-dependence in the color 

^We're concentrating for the moment on the rightmoving nucleus. 

''We are assuming that one is here interested in the field around y = 0, otherwise all the scales in 
the following argument must be boosted appropriately. 



March 9, 2010 ijmpereview 



6 T. Lappi 

charge density and only take the (5-function hmit later at the appropriate stage. In 
the leading logarithmic kinematics the x~ -dependence will be identified with the 
dependence of the color charge density on the rapidity (or x) at which it is probed. 

The computation of the glasma fields from the current ([2]) proceeds in the follow- 
ing wajff^. One must first solve the problem for the field of one individual nucleus; 
the gauge fields of the individual nuclei will then give the initial condition for com- 
puting the field configuration in a two-nucleus collision. We shall first go through 
these small manipulations that will give us a picture of the initial glasma field con- 
figurations after the collision before discussing the relation to the renormalization 
group evolution of the sources. 

The solution for the color current of one nucleus 

J^=5^+pi(xT,a:-) (3) 

is most easily found in the covariant gauge S^A^^^ — 0. One can find a solution 
with only one component of the gauge field is nonzero, namely A'^^^(xt,x~). In 
this case Eq. ^ becomes a 2-dimensional Poisson equation 

- V^A+, = p(xT,x-), (4) 

for which we can formally write the solution as 

A+, = -p(xT,x-)/V^. (5) 

Note that there is an infrared singularity in Eq. ([5]) . The most natural prescription 
to invert the Laplacian is to impose the constraint J d^XT/o(xT) = 0, i.e. to 
require that the source as a whole is color neutral. Imposing color neutrality at a 
shorter length scale will also remove this ambiguity. 

The covariant gauge solution has the advantage of being localized on the light 
cone in the i, z-plane, but its interpretation in terms of partons is not very clear. 
To interpret the classical field in terms of quasi-real Weizsacker- Williams gluons we 
must transform the field into the light cone (LC) gauge. This gauge transformation 
can be done using the path ordered exponential 

U{y:T,x~) ^Pexplig J d?/~A+^(xT, y") L (6) 



givmg 



^^c = (7) 



ALc = -f/(xT,^")a,C/t(xT,x-). (8) 
9 

The light cone gauge solution is not localized on the a; "'"-axis, unlike the one in 
covariant gauge. Instead, it extends to the whole region > as a transverse 
pure gauge field. The field strength tensor F^" , however is nonzero only on the 
light cone x~ = because there is a nonzero energy density only following the color 
source. 
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In general the color charge densities are stochastic random variables drawn from 
some distribution T/Fj,[p]. This probability distribution cannot be computed from 
first principles, but its dependence on y can; it is given by the JIMWLK renor- 
malization group equation, which we shall return to shortly. For many phenomeno- 
logical applications and the following discussion, it will be enough to consider the 
distribution defined by the McLerran-Venugopalan (MV) modeJ^, which could be 
considered as a reasonable initial condition for JIMWLK. In the MV model the 
color charges are distributed with the Gaussian probability distribution 

(p'^(xT, x)p\yT,y)) = g^ti\6'^'5'i^T ~ yT)S{x- ~y-). (9) 

One can motivate the simple Gaussian approximation by an argument that is es- 
sentially based on the central limit theorem. The classical color charge density in 
the CGC is conceptually a sum of the color charges of all the higher x partons 
that have been integrated out of the theory. If these color charges are uncorrelated 
(incoherent) and there are sufficiently many of them, the resulting distribution will 
be GaussiaJS! and local in transverse space. This assumption of independent color 
charges adding up to the effective charge p(xt) also leads to the assumption that the 
parameter should be proportional to the thickness of the target (i.e. ~ A^/'^)] 
a dependence that ismodified by quantum evolution that introduces correlations 
between the chargeJ^^^^. The MV model has two important properties that we 
must mention here. Firstly, as originally shown in it leads to saturation of the 

unintegrated gluon distribution (and of the dipole cross section), at transverse mo- 
mentum scales ^ Qs ~ ofufwe shall comment on the precise relation between 
(7^/1 and Qs later, see also l -^^^ l ^^ i). Secondly, in the dilute large transverse momen- 
tum limit it leads to an unintegrated gluon distribution that behaves as ~ 1/k'^. In 
terms of the integrated distribution this means that xG{x,Q'^) ~ InQ^, the kind 
of generic behavior one expects from DGLAP evolution. These two features; satu- 
ration and a perturbative-like gluon spectrum at large ky, make it reasonable to 
expect that the MV model should be a realistic one for phenomenology at RHIC 
energies, where the effects of JIMWLK evolution have not yet drastically modified 
the distribution of the color charges. 

2.2. Glasma fields 

Let us then turn to the case of two colliding nuclei. It is more convenient here to 
work in an axial gauge. The reason for this is that we would like to understand the 
solution of the classical equations of motion inside the future light cone r > as an 
initial value problem involving only the initial conditions of the fields at r = 
and subsequent propagation into the vacuum. If one studies the solution of the 
two nucleus problem in the covariant gauge the large ^"'"-component of the field of 
the rightmoving source causes a precession (due to covariant current conservation 
[I?^, J^] = 0) of the current of the leftmoving nucleus. Thus in the covariant gauge 
an essential part of the interaction involves the currents directly. In a light cone 
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Fig. 1. Color fields in spacetime. In regions (1) and (2), where only one of the nuclei has passed 
by, the field is the pure gauge field of this one nucleus. In region (3) the field is known numerically. 

gauge one can take advantage of the fact that the solution for the one nucleus case 
has, besides the gauge condition A'^ — 0, also the propert}!! = 0. This is due 
to the static nature of the current; a 2; "'"-dependence in the Wilson line ^ would 
induce a A~ -field. In other words this is a feature of the leading Inx kinematical 
regime, where the degrees of freedom in the source are assumed to have a large 
virtuality ~ p|, ~ Qs^ and therefore a very small p~ . Choosing purely yl+ = 
or A~ = as our gauge condition would, however, break the symmetry between 
the two nuclei in an inconvenient way (see, h owever^ for a discussion of the CYM 
problem in an asymmetric gauge anci^niH] for a formulation of the problem in a 
way that is completely asymmetric from the beginning). Fortunately a convenient 
symmetric axial gauge condition is provided by the temporal, or Fock-Schwinger, 
gauge in the proper time-rapidity coordinates: Aj- — 0. The other independent linear 
combination, A,, — ^t'^A^ = x'^ A^ — A'^ is zero for the one-nucleus solution, 
but not in the case of two colliding nuclei. The Schwinger gauge condition in light 
cone coordinates is x^ A~ + x^ A~ = 0. In other words, on the s"*" = 0-light cone 
where there is a current J^, the gauge condition guarantees that A'^ = and the 
current does not precess, and symmetrically for the other light cone. The effects of 
the interaction between the two nuclei are now visible only in the gauge field, not 
the current (at least at the classical level). 

'^This is for the rightmoving nucleus, for the Icftmovcr = is the light cone gauge condition 
and A+ = a property of the particular solution. 
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Fig. 2. Components of the gauge field, computed numerically on a 512'^-lattice with g^fiR^ = 67.7. 

Inside the future hght cone (the spacetime structure is illustrated in Fig. [T]) the 
gauge fields satisfy the equations of motion in vacuum. What is needed is the initial 
condition for solving these equations. These initial conditions can be obtained by 
requiring that the fields in the different regions match smoothly on the light cone. 
In practice this is done by inserting the ansatz 

A' = e{-x+)e{x-)A\^^ + 0{x+)e{-x-)Ai^^ + 0{x+)e{x-)A\^^ (lo) 
A^^e{x+)e{x~)Al^ (11) 

into the equation of motion ([Ij and requiring that the singular terms arising from 
the derivatives of the 6'-functions cancel. In this way one gets the following initial 
conditions for the gauge field in the future light cone: 

^(3)It=0 = ^(1) + ^(2) (12) 

A;'3)|.=o = |[Aa),A|2)]- (13) 

The equations of motion with these initial conditions can then be solved either 
numerically or perturbatively in the weak field limit. We shall turn to these solutions 
later, but the basic features of the glasma fields can be seen already from the initial 
conditions. 

For a physical interpretation it is useful to calculate the corresponding chromo- 
electric and -magnetic fields in the usual i, z-coordinate system. The field configura- 
tions Eqs. (|12|) and (|13|) correspond to vanishing transverse fields and longitudinal 
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components given by 

^ig[Al,^,Al,^] (14) 

^zge^^[^^,^,Al^]. (15) 

Note that the gauge potential is a true vector and therefore the longitudinal 
component A,, corresponds to Az only at 77 = 0. The field strengths E and B, 
on the other hand are components of the tensor F^"^ and Eqs. ([T4l) and (|15p are 
true for all rj. Solving the equations of motion forward in time then generates also 
transverse components for the fields within a time ^ 1/Qs- A plot of the transverse 
and longitudinal color field strengths as a function of r from a numerical calculation 
is shown in Fig. [2] 

The color fields of the two individual nuclei are transverse electric and magnetic 
fields on the light cone. Why then are the glasma fields in the region between the 
two nuclei longitudinal along the beam axis? One way of understanding these field 
configurations is the following. Let us work still in the Ar = gauge, so that each 
nucleus, when going past a point on the beam axis with no gauge field before the 
collision, leaves behind it a pure gauge field (see Fig.[T]). One can define an effective 
chromoelectric and chromomagnetic charge density by separating the nonlinear 
parts of the vacuum Gauss law and Bianchi identities 

[Di,E']=0 and [D„ B'] = (16) 

as 

d,E' = p, = ig[A\E'] and d,B' = = ig[A\ B% (17) 

Now we can interpret the interaction of the WW chromoelectric and -magnetic 
fields of the nucleus on the x^-light cone with the pure gauge field left behind by 
the other nucleus as an effective chromoelectric and -magnetic charge density left 
behind on the light cone. An exactly opposite charge density is left behind on the 
other sheet, leading to a longitudinal chromoelectric and -magnetic field between 
the sheets. This structure is illustrated in Fig. [3] This description is very similar 
in spirit to the original suggestions for a string based description of a hadronic 
collisions. These string picture had only sources of electric charge ^22]^ When these 
longitudinal field configurations are actually derived from the QCD equations of 
motion also a magnetic field naturally appears. 

It is interesting to note the structure of the energy momentum tensor T^j, = 
jQiivF"^ Pai3 — Ffi^Fva for this initial condition. It is diagonal and, as always in 
gauge theory at the classical level, traceless: T^^ = ^{E^ + B'^) x diag(l, — 1, — 1, 1). 
This can be compared to the standard form for a system with an anisotropy in 
the z-direction: T^i^ = diag(e, — pj_, — p^, — p/,), where e is the energy density and 
p± and PL are the transverse and longitudinal pressures. We see that the initial 
field configuration has negative "longitudinal pressure" . The configuration that is 
the starting point for studies of isotropization by plasma instabilities, where the 
diagonal elements of T^i, at ry = are (e, — e/2, — e/2, 0), is only reached at times 
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Fig. 3. The WW fields of the two nuclei before and after the collision. Before the collision there are 
only transverse fields on the sheets. After the collision the interaction of these fields with the pure 
gauge field of the other nucleus leaves behind an effective electric and magnetic charge density 
(the dots on the figure) on the sheet, and a longitudinal electric and magnetic field between these 
effective charges. 



T > 1 /Qs when the classical fields start to behave linearly due to the expansion of 
the system. 

The glasma field configurations depend on the transverse coordinate at the 
length scale 1/Qs, which corresponds to the typical transverse momentum of the 
gluons being ~ Qs- The same scale is also in general the correlation length of the 
system, and will as such determine the strength of multigluon correlations as we 
shall see in Sec. [51 This leads to the lifetime of the purely longitudinal field con- 
figuration also being ^ 1/Qs- For the bulk of particle production there is therefore 
no clear separation of timescales that would justify treating the field as constant 
in time and space, no matter how appealing this approximation migh be for those 
used to low energy string phenomenology. The difference between the glasma and 
the Lund string model is that conceptually the transverse scale of the problem is 
a semi-hard scale Qs, not the confinement scale Aqcd- Note also that the initial 
fields being longitudinal along the beam axis direction is in no contradiction with 
the lowest order perturbative description of the process as gg ^ g scattering, be- 
cause the longitudinal (with respect to the beam axis) fields are perpendicular to 
the momentum of the gluon being produced. The initial polarization state of this 
gluon is, however, a very particular one. 

T he field i nside the future light cone can then be computed eith er numeri- 
cally I 23 | 24 | 25 | analytically in different approximations (see e.g. for recent 
work). The obtained result is then averaged over the configurations of the sources 

with a distribution Wj,[J'^] that includes the nonperturbative knowledge of the 
large x degrees of freedom. The resulting fields are then decomposed into Fourier 
modes to get the gluon spectrum. This is the method that we will refer to as Clas- 
sical Yang-Mills (CYM) calculations. Note that the average over configurations is 
a cla ssical a verage over a probabilistic distribution. This is guaranteed by a the- 
orem l ^6 | 27 | ensuring the factorization of leading logarithmic corrections to gluon 
production into the quantum evolution of W^y[J''], analogously to the way leading 
logarithms of Q^ are factorized into DGLAP-evolved parton distribution functions. 
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In the hmit when either one or both of the color sources are dilute (the "pp" 
and "pA" cases), the CYM calculation can be done analytically and reduces to a 
factorized form in terms of a convolution of unintegrated parton distributions that 
can include saturation effects: 

diV 11/" d^kr „ , , 

d^PT dy as pf J (27r)^ 

Although this approach (known as "KLN" after the authors o f^, seM for the 
original work) is not strictly valid for the collision of two dense systems, it does have 
the advantage of offering some analytical insight and making it easier to incorporate 
large- a; ingredients into the calculation. We will discuss these calculations of the 
gluon spectrum in more detail in Sec. 21 

2.3. Factorization 

In the previous discussion we did not specify the probability distribution of the 
color charge densities more precisely, besides mentioning the MV model as a phe- 
nomenological approximation. The color charge distribution Wy [p] includes nonper- 
turbative information about the large-x part of the nuclear wavefunction, i.e. the 
valence degrees of freedom boosted from the rest frame of the nucleus. Thus the 
probability distribution cannot be computed from first principles in weak coupling. 
Nevertheless, by considering quantum corrections to high energy scattering, it is 
possible to de rive a renormalization group equation that is known by the acronym 
JIMWLkPSI (pronounced as "gym-walk"). The JIMWLK equation describes the 
dependence of the probability distribution on y, the rapidity (or, equivalently 

in the leading log high energy kinematics that we are working in, on In 1 /a;). In a 
mean-field approximation valid in the large A^c-limit the J IMWL K evolution of the 
correlation function of two Wilson lines reduces to the BkI^SEII equation. 

The role of Wy[p] is analoguous to the conventional parton distribution func- 
tion; it is a nonperturbative quantity whose dependence on one of the kinematical 
variables of the process is described by a weak coupling renormalization group equa- 
tion. In the case of parton distributions, the renormalization group equations are 
the DGLAP ones, and they describe evolution in Q^. In the case of DGLAP one is 
dealing with a dilute system, where the appropriate degrees of freedom are individ- 
ual partons with a definite momentum, whereas in the case of the CGC the good 
degrees of freedom are color charges resulting from interactions of many partons. 
The distributions are similar to parton distributions in the sense that they 

are not (complex) wavefunctions but (at least loosely speaking) real probability 
distributions. This is guaranteed by factorization theorems stating that there is no 
interference between the dynamics of the hadrons or nuclei at large x (JIMWLK; or 
BFKL) or at small (DGLAP) and the process one is studying. Factorization can 
be understood as a statement that one has found the right set of degrees of freedom 
in which one can compute physical observables from only the diagonal elements of 
the density matrix of the incoming nuclei. 
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In our nonhnear high-energy context the JIMWLK equation was derived in 
the context of deep inelastic scattering off a nuclear targe t. In the case of particle 
production in the glasma, it was shown more recently I '^^l^^ l that when one computes 
the NLO quantum corrections to a given observable in the Glasma, all the leading 
logarithmic divergences can be absorbed into the RG evolution of the sources with 
the same Hamiltonian that was derived by considering only the DIS process. The 
underlying physical reason for factorization is that this fluctuation with a large A:"*" 
requires such a long interval in x"*" to radiated that it must be produced well before 
and independently of the interaction with the other (left moving and thus localized 
in x"^) source. 

In the case of the spectrum of gluons produced in a high energy collision the 
leading log divergence can be written as two JIMWLK Hamiltonians, one for each 
nucleus, acting on the expression for the leading order spectrum. This same Hamil- 
tonian describes the RG evolution of the source distributions Wj,[p]. It is most 
naturally expressed as 



■n = \j d'^TA'yTDa{^TK\^T,YT)Dt,{yT) 



(19) 



in terms of Lie derivatives Z?q(xt) operating on the Wilson lines, Eq. ©, con- 
structed from the source color charge densities. The kernel in Eq. (|19p is a function 
of these same Wilson lines: 



»7"''(xT,yT) 



1_ f (xr - Ut) ■ (yt - Ut) 

■K J ^ (xT - UT)2(yT - vltY 



Ut 



U{^T)U\yT) 
~ U{uT)U\yT) 



- ah 

1 . (20) 



The fact that no other dive rgent terms appear is the proof of factorization; this is 
the central result of Ref.^. 



3. Nuclear wavefunction 
3.1. Deep inelastic scattering 

The most direct probe of the nuclear wavefunctions would be to perform nuclear 
deep inelastic scat tering experiments at high energy such as at the Electron-Ion- 
CoUider (EIC) 1331 or a Large Hadron-electron Collider (LHeC) 1311. In the mean- 
time one has to rely on the few existing, relatively low energy, measurements and 
theoretical extrapolations from proton data from HERA. For a review of the ex- 
isting u nder standing of nuclear shadowing we refer the reader to the review by 
Armesto 133. 

It is useful to think of deep inelastic scattering at small x in the dipole picture 
| 37 | 38 | -^jjgj-g process is viewed as a virtual quark fluctuating into a color dipole, 
which then probes the wavefunction of the target (see Fig. |4]). 

In the dipole model one factorizes the total cross section into the probability for 
the virtual photon to fluctuate into a qq pair (a color dipole) and the cross section 
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Fig. 4. In the dipole frame the incoming virtual photon splits into a quark-antiquark dipole of 
trensverse size ry, which then interacts with the target with the dipole cross section a. 




Fig. 5. Left: Wilson line correlator multiplied by fc|, in the MV model, with different discretizations 
of the longitudinal coordinat JE! (see Eq. II25I I below) . Reading off the maximum of these curves 
and defining the corresponding fcy as Qg is one possible way to relate the parameters of the MV 
model to Qs. Right: estimate of the numerical value of the saturation scale, based on a fit to 
HERA data and a Woods-Saxon paramctrization of the nuclear geometrjESE] 

of the dipole scattering with the targel^. The total cross section can be written as 
I4TI581 

aT.i^{x,Q^)^ j d?YT dz|V'T,L(2,rT)|Vdip(x,Q^rT). (21) 

Here the photon wave Junction 'ipT,L{z,rT) gives the probability for the virtual 
photon (T and L stand for, respectively, transverse and longitudinal polarizations 
of the photon) to split into a color dipole of transverse size • The wave function 
'0t,l(^7 tt) includes the known QED part of the reaction and is known analyticalljQ. 



"^There are some tricky issues related to the Lorentz fr ame in which one should view the scattering 

process in the dipole model, see e.g. the discussion inESHSI 

^ To leading order in aom, which is quite sufficient in this context. 
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The exact expressions can be found in e.g. SB. In the classical field approximation 
the dipole cross sec tion can b e expressed in terms of the Wilson lines U{xt) ap- 
pearing in Eq. © ai^OHlESl 

f7dip(rT) = J- / d^bTTr ^1 - U^hr + irT)C/(bT - ^ry)^ 

= 2 1 d^hTAf{hT,rT), (22) 

where we have denoted the imaginary (and dominant at high energy) part of the 
dipole-target scattering amplitude with J\f{hT, tt)- Equation ((22|) provides a direct 
and explicit connection with between DIS observables and the classical field de- 
scription of the initial color fields in a nucleus-nucleus collision. For small dipoles 
the scattering amplitude is proportional to r^, whereas for a large one it approaches 
the unitarity limit J\f ~ 1. The saturation scale Qs is defined as the characteristic 
momentum scale separating these two regimes. Fourier-transforming the Wilson 
line correlator w.r.t. into momentum space Qs could be defined as the cor- 
responding to the maximum of the function k^J\f{hT,^T)- In coordinate space 
one can define it from the value of tt when the scattering amplitude reaches some 
characteristic value. For example one can choose to define QsO^t) by the relatiorP^l 
Af{hT,rT = V2/Qs) - 1 - or h^^Afihr^rT = 1/Qs) = 1 - e^^/^. Several 

other constants in lieu of 1 — e~^/^ or 1 — e^^^^ are used by different authors, lead- 
ing to slightly different numerical values. The three cited above have the advantage 
agreeing with the convention used in the well-known "GBW" fit^-*^ to HERA 
data which is discussed below (see Eq. (|24|) ). The momentum space Wilson line 
correlator in the MV model is shown in Fig. [5] (left). 



3.2. Proton-nucleus collisions 

Another avenue to access properties of the nuclear wavefunction is to study proton- 
nucleus collisions and compare them with nucleus- nucleus ones. Here the idea is 
to treat the proton as dilute probe scattering off the CGC of the nucleus. Among 
CGC theorists it is a common terminology to denote generally the (formal) dilute- 
dense limit of scattering by "-pA^W '^^^^ '^Qj-^ many of the results are equally valid 
for forward scattering in AA collisionJ^. Deuteron-nucleus (falling into the same 
dense-dilute cathegory) collisions at RHIC have turned out to be a very powerful 
tool to access properties of the small x nuclear wavefunction. Indeed the suppression 
of high pt hadron production in pA collisions at forward rapidities observed by the 
BRAHMS collaboratioEp2l have been considered as one of the clearest direct exper- 
imental signals from RHIC favoring gluon saturation. Much of the theoretical con- 
text and the earlier phenomenologial applications have been extensively presented 
in the review by Jalilian-Marian and Kovchego\E2l, and we shall here discuss them 
only briefly. More recently the focus in pA collisions as well as nucleus-nucleus ones 
has been on azimuthal and rapidity correlations; we shall return to these towards 



March 9, 2010 ijmpereview 



16 T. Lappi 



a 
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Fig. 6. Ha rd particle spectra in dAu collisions using the DHJ parametrization of the dipole cross 
sectionJ3S]. 



the end of this section. 

In pA-colhsions the particularly interesting kinematical regime is in the forward 
rapidity region (proton fragmentation region, i.e. small x in the nucleus and large x 
in the proton) , where the nuclear saturation scale is large and the proton (deuteron 
in practice) is a dilute probe. Here one is typically dealing with particle production 
at large transverse momenta compared to the intrinsic in the proton wavefunc- 
tion, and it is necessary to resum the large DGLAP logarithms on the proton side. 
One must also take into account both the quark and gluon degrees of freedom from 
the proton, whose scattering off the CGC target depends on the Wilson line cor- 
relator in different representations. The single inclusive hadron spectrum can be 



written aa 



dcrP^- 



■hX 



d2/d2pTd2bT (27r)2 



da;. 



/i 

fg/p 



Tq/p\Xp 



ixp,Q^)M^ 



XF 



-Pr, 



D 



— Pt,ot 

Xf 



D 



h/g 



h/q 

Xp 
Xr> 



Xp 
Xp 



} (23) 

Here fq^g/p are the gluon distribution functions in the proton and Dh/q,g the frag- 
mentation functions of quarks and gluons into the hadron h. The Wilson line cor- 
relators Af,a (as in Eq. ([22| ) in the fundamental and adjoint representations have 
to be evaluated in at the rapidity scale {x the nuclear wavefunction corresponding 
to the produced hadron. A similar approximation for only the gluonic contribution 
can be derived from the /cr-factorized formulation discussed in Sec. |4] and has also 
often been used in applications to the spectrum in pA collisions. A comparison 
of Eq. ([23|) using a certain parametrization ( "DH J"I4£I) with experimental data is 
shown in Fig. |6l Computing the spectrum in pA collisions at different rapidities 
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r^^^ 
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Fig. 7. Left: nuclear modification factor R^ a fro m BRAHM352. Right: quaUtative prediction of 
the evolution of R-pA from a CGC calculatiorlSfii 



enables one to almost directly see the effects of high energy evolution in the nu- 
clear wavefunction on the Wilson line correlators A/f^a- At the starting scale of 
the evolution one generically expects a "Cronin" enhancement in pA compared to 
pp collisions. This, as can explicitly be seen e.g. in the M.V. model, is due to the 
saturation in the nuclear wavefunction that suppresses the gluon distribution below 
and causes an enhancement above (the "Cronin peak"). The effect of high 
energy evolution is to wash this peak awaj^^^"^^-^, which is seen in the data at 
more forward rapidities, as shown in Fig. [71 

A more r ecent area of activity are two-particle correlations in dAu 
collisionJSlEi. The baseline comparison is the structure observed in proton-proton 
collisions, where high particles typically originate from back-to-back jets. The 
azimuthal structure of the two particle correlation is thus a rather narrow peak 
on the away-side resulting from the fragmentation of the original jet into a narrow 
cone. The expectation from the CGC picture is essentially that of an increased 
collectivity in the cor relation, signaling itself as a flattening and broadening of the 
away-side peali™. This feature should become more prominent as the kinemat- 
ics of the two particles move towards small xa (when the saturation scale becomes 
larger) and should be very sensitive t o th e relation between the trigger pt's and 
Qs- A result of a particular calculatiorf^ is shown in Fig. [51 together from recent 
experimental data from STAlP^. 

3.3. Dipole cross section parametrizations 

By now there is a wide variety of saturation-based parametrizations of the dipole 
cross section, mostly fitted to HERA data. While the aim here is not to give 
a comprehensive listing, let us mention a few. One of the problems of practical 
parametrizations is that calculations of DIS observables (as in Eq. (|2T|l ) use the 
dipole cross section as a function of ry, the size of the dipole, what one needs 
to compute particle production in pA collisions (as in Eq. (|23|)) is the Fourier- 
transform. Although it is not strictly speaking formally required, one would prefer 
to have a parametrization that is positive definite in both position and momentum 
space, while retaining some other general features required in d iffer ent limits. It 
has turned out to be a surprisingly difficult mathematical problerrP2l to find such a 
parametrization. The result of difficulty has been that many authors prefer to use 
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Fig. 8. Predictio n fo r broadening of the two-particle correlation function in deuteron-nucleus 
collisions from Ref.ED. Right: correspondinf data from the STAR collaboration^^. 



different parametrizations for DIS and for pA studies, thus reducing the predictive 
power of the respective calculations. 

The early and widely used "GBW'SlSSl dipole cross section is a Gaussian in 
yt, and consequently also in kj^: 



o'dip(a;, Q^jTt) = cro(l - e 



)/4)^ 



(24) 



with Qs^{x) ^ with A ~ 0.3 resulting from the fit to HERA data. This for- 
mula features the idea of parton saturation very clearly. The exponential decrease 
at high momenta is, however, not very physical. At very large k-p one would ex- 
pect CTdiplkr) to behave as ~ l/kf^ in order to recover a DGLAP-like increase of 
the gluon distribution at high momenta: xG{x,Q^) ^ InQ^; one of the attrac- 
tive features of the MV model is precisely that it results in this behavior. At 
somewhat smaller momenta, kx ^ Qs BK/JIMWLK evolution predicts a differ- 



ent power law crdip(kT) ~ k, 



-2(1+7) 



with an anomalous dimension 7 0.6. In the 



BKW" parametrizatiorP2l the large kx behavior is modified by replacing the 
ponent r^Qs^ by r'^xG{x,l/r^), where xG{x,l/r'^) is a DGLAP-evolved gluon 
distribution evaluated at the scale /Iq -I- 4/r|,. This brings the high-fc^ limit closer 
to a DGLAP-like behavior, but applying the DGLAP-evolved distribution around 
the saturation scale is not without problems; in particular it turns out that the best 
fit to data is given by a gluon distribution which decreases towards smaller x at 
the initial scale A dipole cross section that reproduces, instead of DGLAP, the 
behavior of BK evolution in the extended scaling region kx ^ Qs is given by the 
"IIM"IS1 parametrization. 

Another simplification in Eq. (j24p is the dependence on the impact parame- 
ter hx- The impact parameter profile of the dipole scattering amplitude can be 
directly measured from the i-distribution of diffractive events; experiments are 
consistent with a Gaussian in hx (implying an exponential dependence on t). In 
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Eq. the impact parameter dependence T(bT) has been factorized from the yt- 
dependence, which reduced it in the total cross section into a constant proton area 
ctq = / d^hxTiyiT)- The problematic aspect of this factorization is that the unitar- 
ity limit of the sc att ering amplitude is only approached at 6 = (see more detailed 
discussion in RefP^. Two parametrizations attempting to treat the 6-dependence 
more consistently, either in a DGLAP-improved approach or parametrizing features 
of BK evolution are provided by the IPsat (also "KT'^ ) and bCGC models^?'. In 
these parametrizations the saturation scale is impact parameter dependent. 

The "KKT"ESland "DHJ"HS1 parametrizations, also including a BK-like anoma- 
lous dimension, have been used in pA collisio ns, but not fitted to DIS data. There 
are also by now several numerical evaluation j^'^ l ^^ of the BK equations numeri- 
cally both with fixed and running coupling. The initial condition of the evolution 
can then be used to fit experimental data and produce a parametrizatioiJ^ for use 
in other contexts. A recent application of a direct solution of the BK equation to 
hadron spectra in pA collisions can be found in Ref.-A 



4. Bulk gluon production in nucleus-nucleus collisions 

Let us now turn back to the case of the collision of two dense systems. The first 
question to understand are the properties, such as entropy or energy density, of the 
initial stage in a heavy ion collision in the CGC framework. As we have discussed, 
the initial stage of the collision system at central rapidities is, to a first approxi- 
mation, gluonic and boost invariant. It is therefore far from chemical and thermal 
equilibrium (in particular anisotropic, with much larger transverse than longitu- 
dinal momenta) . The thermalization of this gluonic system towards a quark-gluon 
plasma is not quantitatively understood, in spite of a lot of recent work on exploring 
the instabilities that could dominate this stagJ^. 

To make a more direct connection to phenomenology one must thus at this stage 
try to shortcut the thermalization stage by replacing it with a more schematic 
relation between the initial gluon spectrum and the later stage evolution of the 
system. One simple version of this relation, based on experimental observation, 
goes by the name of "parton-hadron duality" , stating that the final state hadron 
multiplicity is proportional to that of the initial partons. A second way of estimating 
the relation is based on the argument that the particle multiplicity cannot decrease 
during the thermalization process, due to the second law of thermodynamics and the 
fact that the entropy is essentially proportional to the particle density. On the other 
hand the energy per unit of rapidity can only decrease (due to pAV work, i.e. fiow 
of energy towards the fragmentation region) during the time evolution. Thus one 
can to some extent estimate the relation between the initial gluon multiplicity and 
energy density and the final (charged hadron) multiplicity and transverse energy. In 
particular it seems reasonable to assume that the dependences of both on centrality, 
collision energy and rapidity would be similar. 
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4.1. Calculating the initial gluon multiplicity 

It is thus of direct phenomenological interest to try to calculate the initial gluon 
multiplicity in the CGC framework with, as much as possible, input values of the 
parameters determined from independent observables, in particular from DIS or pA 
data. As already mentioned in Sec. [21 there are two major ways of doing this. One is 
to directly solve the classical Yang-Mills equations numerically ( "GYM" , the other 
is to use a fcx-factorized approximation ("KLN") in a situation where it strictly 
speaking is not valid. The results from both approaches for the integrated gluon 
multiplicity, if not for the spectrum as a function of momentum, turn out to be very 
similar. This can mostly be understood by dimensional analysis. Gluon production 
in a central collision, at midrapidity, is a one scale problem with Qs as the only 
relevant dimensionful scale. The gluon density must therefore be dN/ dy d^x^ ~ 
Qs^. In a noncentral collision or away from midrapidity the saturation scales of 
the two nuclei are different, but the generic result in both the GYM and KLN 
calculations is that dN/ dy d^XT ~ Qsmin' where Qsmin the smaller of the two 
saturation scales. Thus the dependence of the integrated multiplicity on centrality 
is determined by the impact parameter dependence of the saturation scale, and the 
dependence on rapidity and collision energy by its x-dependence. The next natural 
observable to look at would be the energy density, which probes an integral of the 
gluon spectrum weighted by an extra power of px- This is more sensitive to the 
large momentum part of the spectrum, which would typically behave as 1/p^ (as 
resulting from the expected l/fc|,-behavior of the unintegrated gluon distribution) 
in both cases. The actual gluon spectrum from the GYM and KLN calculations can 
be very different especially for kx ^ Qs, but it is difficult to get an experimental 
handle on the shape of the spectrum. 

Typically in KLN calculations, if the value of Qs is fixed by an external input, 
the normalization of the result is adjusted by hand to a reference point in the 
data. In particular, in the A:T~factorized approximation in the cases where one can 
derive it from the GGG framework the integrated gluon spectrum is not IR finite; 
although the divergence is only logarithmic instead of a power law due to saturation. 
In practical calculations one must thus add an additional regulating prescription 
with the corresponding adjustment to the normalization. This adjustment of the 
normalization is not possible in the GYM calculation, where the only parameter one 
has is essentially Qs, which determines both the number and the typical transverse 
momentum of the gluons produced. Thus in the GYM framework it has turned 
out to be more realistic to make a genuine prediction for the normalization of the 
spectrum based only on the value of Qs obtained from HERA data; we will discuss 
this estimate below. Due to the relative simplicity of the analytic fcr-factorized 
approximation, on the other hand, it is easier to incorporate more detailed features 
of the transverse coordinate or rapidity dependence. 

Most of the original GYM calculations were done using the MV model, treating 
the color charge density g^p as a free parameter to be adjusted to RHIG data. We 



March 9, 2010 ijmpereview 



Small X physics and RHIC data 21 

shall here follow the approach of Refff^ to take these results and transform them a 
posteriori to parameter-free RHIC postdictions using the numerically determined 
relation between g^fx and Qs and values of Qs extracted from HERA fits. 

Let us first recall the results of the numerical CYM calculations of the gluon 
spectrum in heavy ion collisions ^"^ '^'^ '^^ We shall not discuss the numer- 
ical procedure here; a recent review is given in^^. In order to compare the results 
we must, however, comment on one aspect of the numerical implementation of the 
MV model. In the numerical implementation of the MV model the Wilson lines ^ 
are constructed as 

f 1,2/ N ^ 

U{.t) = n -'^-"^^tH ' (25) 

1 1 I + ™' J 

where the color charges are Gaussian variables with the variance 

(p^(xT)pf(yT)) - S^^'S'^'S^i^T - yr)^. (26) 

The indices k,l — 1 . . . Ny represent a discretization of the longitudinal direction 
into Ny small steps; the continuum limit corresponding to Eq. ([6]) is achieved for 
Ny — >■ cx) at constant g^i^A- Some kind of infrared regulator is needed in order to 
invert the Laplacian operator V ^, it can be regularized by the scale m in Eq. ([25l) . 
As pointed out in Ref.l^(see alscH^, the relation between the parameter g^fj, of the 
model and the physical length scale Qs, the correlation length of the Wilson lines 
in the transverse plane, depends on Ny. The conceptual picture CGC framework is 
that of a Wilson line built up from infinitesimal steps in rapidity, i.e. Ny — >■ oo. The 
early numerical implementations, on the other hand, used Ny — 1 for simplicity. 
When expressed in terms of the physical length scale Qs the results depend very 
little on the way the Wilson lines were constructed. However, when comparing 
the results of the numerical calculations to phenomenology and using values of Qs 
extracted independently one must take care to use the correct relation between Qs 
and g^fj,. As shown irflSi, the relation corresponding to m — >■ 0,Ny = 1; i.e. the 
prescription used in the numerical CYM calculations, is Qs ~ 0.57g^/x. 

With this conversion between Qs and g^^ explicitly stated let us then proceed 
to the result of the numerical CYM computations. The energy and multiplicity per 
unit rapidity can be parametrized as 

^ = 7 ^^^^ 

- -2 (28) 

The numerical result (see in particular Refs. IMUS) jg ~ 0.25 and /^r w 0.3. 
We can also recast this e xpression in terms of the the "liberation coefficient" c, 
introduced by A. Mueller l ^5 | 76 | 77 |^ rpj^^ liberation coefficient is defined by writing 
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the produced gluon muhiphcity as 



dN CfQ 



2 



With Eq. ([271) this leads to 



(29) 



c=g^(9±\ . (30) 



2Cf V Q: 

The original expectation was that c should be of order unitjE^t^. The analytical 
calculation by Y. Kovchegov 1^ gave the estimate c w 2 In 2 « 1.4. With /at = 0.2 
in Eq. ([27| and Qs/ff^M ~ 0.57 we see that the CYM result for the liberation coeffi- 
cient is c sa 1.1. A demonstration of the independence of this number on the details 
of the Wilson line correlator can be made by replacing the MV model distribution 
by another one. For examp le, if the Wilson lines are co nstruct ed explicitly to cor- 
respond to the iPsatl^fiilir^Xowalski-Teaney) or bCGC (EHIl}) parametrization of 
HERA data, the result for the liberation coefficient c is the same as in the MV 
model across the whole energy range from RHIC to the LHCp^l. This is the result 
in spite of the very different forms of the gluon spectra at high fc^ in the models. 

The estimates for the values of Qs based on HERA data, n uclear geometry 
and attempts to fit the limited available nuclear DIS 

dati™ 

vary somewhat, 

especially due to the degeneracy between the proton size and the value of Qs in 
the HERA fits; see e.g.^l^l for a comparison of the numerical values. We take here 
the IPsat estimat^^ that removes this interplay by fixing the size of the proton 
from diffractive data, where it can be directly measured. The estirnateamounts to 
Qs ~ 1.2 GeV for an average central RHIC collision at midrapidit^ ^ 1 6 | 1 8 j noted 
previously, this corresponds to the MV model parameter g^/x w 2.1 GeV (at Ny = 1) 
in the CYM simulations. Plugging this number into Eq. (|27p leads to the estimate of 
ss 1100 gluons in the initial stage of a central gold-gold collision at midrapidity 
at RHIC. This is very close to the observed final total (charged-l- neutral) parti- 
cle multiplicity. In other words, recalling our discussion earlier in this section, this 
points towards a very rapid thermalization that practically conserves particle num- 
ber, and leaves very little room for higher order effects to increase the multiplicity 
during the thermalization stage. 



4.2. Energy and rapidity dependence 

As we discussed, the gluon multiplicity is, across different parametrizations to a 
very good approximation proportional to itR\Qs^ / as- Thus the predictions for 
LHC collisions depend mostly on the energy dependence of Qs- On this front there 
is perhaps more uncertainty than is generally acknowledged, the estimates for A = 
dlnQs^/ dlnl/a; varying between A = 0.29 1^ and A = 0.18 ^'^ in fixed coupling 
fits to HERA data, with a running coupling solution of the BK equation giving 
something in between these values'^. This dominates the uncertainty in predictions 
for the LHC multiplicity (see Fig. [9|). 
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Fig. 9. Predictions for the evolution speed from the running cou pling BK equation and the 
subsequent gluon multiplicity extrapolated to LHC energies from Ref!S2l. 

The RHIC coUision energy is still too slow to clearly see any saturation effects in 
the rapidity dependence of the multiplicity around y — Q. K simple estimate for the 
effects of large x physics, such as momentum conservation, is to consider the typical 
(1 — x)*-dependence of gluon distributions at large x. Inserting x = {p±) / y/s 
leads to the estimate Ay ~ y^S^/s/ {p±) for the rapidity scale at which the large x 
effects contribute to the rapidity distribution around y = 0. This means that the 
(1 — a;)^-behavior starts to dominate the rapidity dependence of the multiplicity 
around midrapidity at a scale of Ay ~ 4 RHIC and Ay ^ 19 at the LHC. Note 
also that the large x contribution is an effect of order 1 at this scale, whereas small 
X evolution can be expected to give a much smaller effect at a rapidity scale 
Ay ~ l/cKs ^ 3. Only at the LHC the large x effects will be mostly absent around 
midrapidity and one has a good possibility of seeing a clear signal of CGC effects 
in the rapidity dependence of the multiplicity. 

5. Transverse geometry 

Relativistic heavy ion collisions can take place with any impact parameter, from 
very peripheral ones that should look like simple nucleon-nucleon collisions to cen- 
tral ones that actually produce a system of the size of the colliding nuclei. Ex- 
perimentally this provides a tool to study quantities as a function of the system 
size, assuming that one is able to detect the centrality. The impact parameter not 
being a measurable quantity per se, there are basically two methods used to ex- 
tract the collision centrality in experiments. One is to measure the noninteracting 
spectator nucleons in the zero-degree calorimeters (ZDC) and the other to make 
a simple assumption on the dependence of the charged multiplicity in some part 
of the detector on the size of the system. This is usually done via a Monte Carlo 
Glauber calculation, which enables one to compute the distributions of A^part and 
-^coii (the numbers of participant nucleons and binary nucleon-nucleon collisions) 
for a fixed impact parameter; assuming that the individual nucleon-nucleon colli- 
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Fig. 1 0. Centralitv dependence of the multiplicity. Left: CYM calculation • ^ Right: KLN calcu- 
lation EH 

sions are independent of each other. Using a simple ansatz for the dependence of 
the charged particle multiplicity on A'part and A'coii one can then divide the events 
of the whole minimum bias data set into centrality clas ses and estimate the typical 
impact parameters corresponding to each class (see for a review of Monte Carlo 
Glauber modeling). 

Besides the ZDC data there are few ways to independently check the consistency 
of the MC Glauber framework for understanding the collision geometry. One is is 
nevertheless comforted by a general impression that the picture seems to work very 
consistently, depending weakly on the details of the MC Glauber model and success- 
fully parametrizing a wide range of phenomena. In particular it turns out that the 
bulk particle production mechanism is such that the charged particle multiplicity 
is roughly proportional to the number of participant ("wounded") nucleons, with 
a constant of proportionality increasing slightly as one goes to more central colli- 
sions. A MC Glauber computation is, however, merely a parametrization of data: 
it does not contain any dynamics and therefore does not provide a microscopical 
explanation of the particle production mechanism. 

5.1. Centrality dependence of multiplicity 

The approximate proportionality of dN/ drj to A^part is a natural result in the CGC 
framework. The saturation scale Qs being the only dimensionful parameter in the 
problem the multiplicity per unit area should be ~ In a fully central collision 

the saturation scales of the two nuclei at a transverse coordinate xt are equal to 
each other and proportional to the nuclear thickness (and therefore to iVpart) at the 
point XT', integrating over the transverse plane gives a multiplicity proportional to 

-Apart • 

For noncentral collisions the situation at a fixed transverse coordinate is 
asymmetric since the two saturation scales and are not equal to each other. 
In this case the multiplicity can also have a dependence on Q^/Q^ . Parametrically, 
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the spectrum of produced gluons in this case behaves l ^^ l ^^ l as 



d^XT d^PT 



In(pT), PT<Qsi (31) 



2 



2/ 



Qsl<PT <Q^2 (32) 



4 , PT > Qs2, (33) 

where Qsi and (5s2 are the smaller and larger of the saturation scales Qf and . 
Integrated over transverse momenta, this gives 

diV o 

d^ - Qs? (34) 

neglecting logarithmic corrections ^ In (Qs2/Qsi)- 

Both the CYM and KLN types of calculationj2Hl84l86| 

reproduce quite well 

the experimental data on the centrality dependence of the multiplicity, as shown 
in Fig. [TOl In the CYM calculation all the deviations from a strict proportional- 
ity dN/ Ar] ~ iVpart come from the (parametrically logarithmic) deviations from 
the strict proportionality to the smaller one of the saturation scales, Eq. (p4| . 
In the KLN calculations, especially in the earlier versions where the geometry 
was treated in a more simplified way, a significant part of the deviation from 
d.N I drj A^part was argued to result from the running of the coupling constant 
in dN/ drj ~ nR^Qs^ /as{Qs)- The deviations from a number of participant scaling 
in the experimental data are small enough that they do not allow one to discrimi- 
nate between the relative importance of these effects. 



5.2. Eccentricity of the initial state 

A striking signal of collective behavior of the matter produced at RHIC is ellip- 
tic flow (for a recent review, seJ^^. Comparing hydrodynamical calculations with 
flow is a way to address fundamental properties of the medium, such as the vis- 
cosity, but this comparison requires understanding of the initial conditions of the 
hydrodynamical evolution, particularly the initial eccentricity for elliptic flow. The 
original general consensus some years ago was that ideal hydrodynamics is in good 
agreement with the experimental data. Among other caveats this general conclu- 
sion also supposes that one is relatively free to choose the initial conditions — in 
particular the transverse geometry — of the hydrodynamical calculation to fit the 
experimental data (see e.g. the thorough comparison of initial conditionJSHl) ^ This 
claim has had to be reevaluated more recently after it was argued (using a KLN-type 
calculation^^ 90 91) ^i^^^ qq(j results in a larger initial eccentricity than tradition- 
ally used in hydrodynamical calculation (mostly participant or wounded nucleon, 
often called "Glauber" initial conditions) . This leaves more room for viscosity in the 
hydrodynamical evolution. The CGC estimate has now settled to a lower value than 
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Fig. 11. Left: The original KLN eccentricity ESI compared to the CYM resultEH. Right: A subse- 
quent revision of the KLN model was to reformulate the nonuniversal saturation sc ale in terms of 
fluctuations in the participant nucleons led to a lower eccentricity estimate ("fKLN")E3practically 
in agreement with the CYM result and very close to binary collisions scaling. 



first argued irls^due to further developments that we will outline below. The result 
remains, however, that generically CGC calculations tend to predict an initial state 
with an eccentricity close to or slightly above iVcou scaling, which is larger than 
often used in hydrodynamics. The eccentricity obtained in the CYM calculation is 
demonstrated in Fig. 1111 Everything else staying equal, the larger eccentricity of 
the initial state can be compensated by the effect of a larger viscosity in order to 
produce the observed elliptic flow V2- When other parameters in the fit are also 
allowed to change, the situation can become more complicated and, in some ca ses, 
contrary to this basic intuition. For example in the "Knudsen number" fit of^ 
assuming a finite cross section (and thus deviations from ideal hydrodynamics in- 
cluding viscosity) also the equation of state was allowed to vary as a fit parameter. 
The result of the fit was that the CGC initial conditions actually corresponded to 
a smaller viscosity (larger cross section) ; the observed smaller V2 being achieved by 
a softer equation of state (cs = 0.22 vs 0.3 for the Glauber initial conditions). 

As we saw previously, in the case of two different saturation scales at a point 
in the transverse plane, the gluon multiplicity only depends on the smaller one of 
the two. The initial eccentricity, on the other hand, is computed from the energy 
densitj0, which behaves parametrically as QsiQs2j where again Qsi < Qs2- the 
original KLN calculations the transverse coordinate dependence of the saturation 
scale was taken as (xt) ^ "■part,A(xT), where npart,A(xT) is the density of 
participating nucleons in nucleus A. This has the advantage of being directly con- 
nected to the well established Glauber modeling of the geometry also used to bin 

' In particular in the CYM calculation, where one is explicitly faced with the fact that the energy 
density is practically the only local gauge invariant observable available. The gluon spectrum is 
determined by Fourier-transforming the field modes and is therefore necessarily nonlocal at a scale 
1/pt- In a (semi) analytical calculation such as KLN one can formally evade this problem. 
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the experimental data. To understand the effect of this prescription on the cal- 
culated eccentricity one must study what happens in a noncentral collision. The 
important region for this is the one at the edge of one of the nuclei (nucleus A) and 
the center of the other one (B). In this case the definition leads to both saturation 
scales being small (since npart,B goes to zero outside of nucleus A). This is concep- 
tually problematic, since one would expect the saturation scale to be large in the 
center of nucleus B. The KLN prescription is therefore nonuniversal: the value of 
Qf is not a property of nucleus B alone, but is determined by a final state effect; 
namely the presence or not of the other nucleus A at a given transverse coordi- 
nate. Since e ^ QfQf in this region, this prescription leads to a suppression in 
the energy density in the edge of the interaction region and a larger eccentricity, 
basically independently of the precise form of the unintegrated gluons distribution 
(see Fig. I12p . To remedy this nonniversality problem in the KLN model it has sub- 
sequently been reformulated into what is known as the fKLpJ^SI (fluctuating KLN) 
model, where one first introduces explicitly the nucleon-nucleon cross section to 
determine whether a given nucleon participates in the scattering. After this step 
the resulting gluon spectrum is the calculated with universal saturation scales in 
the two nuclei. The fKLN model, has practically by construction, the property that 
a) the saturation scales are universal, b) the original KLN centrality dependence 
of the multiplicity is reproduced (due to the property ([34])) and c) it reproduces 
event-by-event fluctuationiSHl that are close to ones given by Monte Carlo Glauber 
calculations. (For experimental results on V2 fluctuations see e.g. Ref.^^Sl), xhe re- 
sulting eccentricity, as shown in Fig. [TT] (right) becomes smaller, closer to the result 
(similar to A^coii scaling) obtained in the straightforward and universal CYM calcu- 
lation. The price to pay for this modification is a nonuniversality of another kind: 
namely the fcT-factorized formalism becomes different for AA-collisions than it is for 
the case of proton-proton collisions. In the treatment of AA collisions in the fKLN 
formulation there is an additional step of computing first a collision probability 
of individual constituents of the colliding objects (nucleons in a nucleus) which is 
not invoked in the pp-case. In other words, nuclei are, even at high energy and 
high pT^ not treated as consisting of quarks and gluons, but as fundamentally more 
complicated objects to which the usual (coUinear) QCD factorization theorems do 
not apply. One curious consequence of this feature is that the saturation scale of 
the nucleus in the fKLN model approaches a constant, not zero, at arbitrarily large 
distances outside the nucleus. 

6. Correlations in the glasma 

So far we have discussed how the CGC picture of the small x nuclear wavefunc- 
tion can be used to compute bulk properties of the initial state of a heavy ion 
collision. Experiments do not probe this initial state directly, because the system 
goes through a complicated time evolution before the hadronization stage. A good 
candidate for a more direct experimental probe of the glasma stage is provided by 
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Fig. 12. Left: energy density as a function of the transverse coordinate in the original KLN 
definition compared to participant number scaling. Right: the KLN eccentricity for different un- 
intcgrated gluon distribution parametrizations. Plots from Ref.l2QI. 




Fig. 13. Left: forward- bac kward multiplicity correlations from STAll241 Right: Experimental 2- 
particle correlation dat jSS] in AuAu and dAu collisions, showing the elongated ridge structure in 
the former case. 



different kinds of correlation measurements, and they have recently been a focus 
of both experimental and theoretical activity. Particles that are produced far in 
rapidity can, by causality, only be correlated at early times, and should therefore 
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be directly related to the properties of the glasm£@. Examples of these long range 
correlations are the "R idge" structure of two particle correlations in central heavy 
ion collisions I96|97|95 | g^^^ lovLg range rapidity correlations in particle multiplici- 
ties fS^'^^i (see Fig. [131 ^- Anothe r interest ing correlation measurement are electric 



charge-reaction plane correlationJl^QlI^ that can be due to CP-violating fluctua- 
tions in the initial stage. As we shall argue, these correlations also arise naturally 
in the glasma. Correlations between the electric charge and the reaction plane are 
most likely mediated by the magnetic field caused by the noncentral collision of 
two positively charged ions; it is hard to think of another quantity that would as 
directly couple to both the reaction plane direction and the electric charge. Since 
this magnetic field dies away very fast with r, the QCD part if the explanation must 
be sensitive to the gluon fields at the earliest stage of the collision. In the following 
we shall discuss some features of how these correlations can be understood in the 
glasma framework. 

6.1. Multigluon correlations in the boost invariant case 



In Sec. l2.3l we argued that the leading logarithmic corrections to particle production 
should be factorizable into the distributions describing the sources. As the argument 
is essentially one of c aus ality, it applies also to multigluon production (a more formal 
proof is given in Ref.l^^ . The case is particularly straightforward when the rapidity 
differences between the produced gluons are small enough compared to l/ckg, so 
that at leading order one does not need to resum the radiation of additional gluons 
between the measured ones. In this case the rapidity structure is indeed a trivial 
one: when observing gluon fields within a rapidity window smaller than I fas', all 
the gluons are sensitive to the same configuration of sources and the gluon fields 
are boost invariant. 

In fact, gluon correlations in AA-collisionf0 are in a sense simpler than in more 
dilute systems. To understand why let us first formally evaluate the leading con- 
tribution in a dense-dense collision. Consider the probability distribution of the 
number of gluons produced in a small rapidity interval. It is convenient to define a 
generating functional 



^ " (35) 



d^pi • • • d^pn ■ 



Y[ d^p, (z(p,;) - 1) 

.2=1 

The Taylor coefficients of J- around z = I correspond to the moments of the 
probability distribution; integrated over the momenta of the produced gluons they 

^Although when making precise estimates one must be more careful and distinguish spacetime, 
momentum space and pseudorapidity and understand shorter range correlations such as those 
from hadronization. 

^ The interpretation of this multiplicity correlation data is somewhat puzzl ing because of the 



interplay with the impact parameter induced, purely geometrical, correlations I2S 
' By "AA" what is meant here is really the formal power counting situation when the color sources 
of both nuclei are strong, as we shall discuss shortly. 
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Fig. 14. Relative importance of connected and disconnected diagrams to the two gluon correlation 
function. One of the color charge densities is considered large, pi ~ whereas the other is 

allowed to vary between the "AA" case p2 ~ l/g and the "pA" one p2 ~ g. The order of the 
disconnected diagram, on top, is g'^p'^p'^, whereas the interference diagram in the middle is g'^p'^p^ 
and the connected one, lowest, is g^p^p^. In the "AA" case the disconnected diagram dominates, 
for the "pA" case all three are equally important. In the dilute "pp" limit only the connected 
diagram matters and both gluons are produced from the same BFKL ladder. 



are 

{N) {N{N-l)) ... {N{N -1)---{N ~n + l)) . (36) 
97 

The result of Ref. ' is that when these moments are calculated to NLO accuracy, 
the leading logarithms can be resummed into the JIMWLK evolution of the sources; 
completely analogously to the single inclusive gluon distribution. The resulting 
probability distribution can be written as: 

VKy pi WV P2 -7-?^ -Ti — e ■> (37) 



d^Pi • • • d^Pn J ^ nl d^pi d^p. 

Pl>P2 

Here the factors denote the multiplicities corresponding to the leading order 

classical fields. 

The probability distribution (157)) has been derived to the leading order in a weak 
coupling expansion in as, and resumming leading logarithms of l/x. In this power 
counting, each insertion of the external color source contributes a factor p ^ 1/5, so 
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that the saturation scale g^{pp) ~ Qs^ is taken to have no powers of the couphng. 
This leads to the glasma fields being parametrically ^ 1 / <? and ^ 1 at 
the dominant classical level. An additional insertion of a color source then does not 
change the power of as of the result, since the l/g in p compensates the additional g 
in the three gluon vertex that attaches it into the diagram. Thus the leading contri- 
bution is given by a s um o f all tree diagrams with an arbitrary number of insertions 
of the classical sourc a computation that has to be done numerically. In the 
language of BFKL physics, the dominant contribution to n-gluon production in this 
"AA" case comes from cutting one rung in n distinct BFKL ladders. Producing two 
gluons from the same ladder is suppressed by the weak coupling compared to this 
"disconnected" contribution. In the "pA" or "pp" power counting, on the other 
hand, one or both of the sources p are assumed to be weak. For example at large x 
the color source consists of a number 0(1) of charged (valence-like) partons, each 
with a charge ~ g (times a color factor). This leads to p ^ g, and additional inser- 
tions of such a weak external color charge lead to a contribution that is suppressed 
at weak coupling. Thus in the "pp" case the dominant contribution to multigluon 
production is achieved by producing all the gluons from the same ladder, mini- 
mizing the number of insertions of the color source. Whereas in the AA-case the 
dominant contribution comes from (complicated) tree diagrams and has a classical 
field interpretation, in the "pp" case the n gluons are produced from the same lad- 
der, which is an n — 1-loop diagram. In the "pA" case, when pA ^ ^/g and Pp ^ g 
both the connected (loop) and disconnected digrams are parametrically equally 
important. This structure is illustrated in Fig. 1141 where the contributions of the 
"disconnected" and "connected" diagrams and their interference term are evaluated 
for the case of pi ~ 1/g and p2 having different parametric dependences on g. In 
practice, however, many applications to two particle corre lations in the "pA" case 
only the connected part has been considered (see e.g. for a discussion). 

Note that the Poissonian-looking form of the result is to some extent an artifact 
of our choosing to develop and truncate precisely the moments Eq. (|35|) that are 
simply (N)"' for a Poissonian distribution. Since in our power counting N ^ 1/as, 
any contributions that would make the distribution Eq. (|37p deviate from the func- 
tional form are of higher order in the weak coupling expansion of the moments (|35|) 
and are neglected in the calculation unless they are enhanced by large logarithms 
of X. It should also be emphasized that in spite of appearances of Eq. ([37|) the prob- 
ability distribution is in fact not Poissonian. To understand the nontrivial nature 
of this result it must be remembered that the individual factors of in Eq. (|37l) 
are all functionals of the same color charge densities pi, 2; thus the averaging over 
the p's induces a correlation between them. These correlations are precisely the 
leading logarithmic modifications to the probability distribution; they have been 
resummed into the distributions Wy] the functional form of the multigluon corre- 
lation function under the functional integral in Eq. ([37]) is the same as at leading 
order. The fact that the dominant contribution is disconnected for fixed sources 
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Fig. 15. Left; Building block, Lipatov vertex coupled to two sources. Right: combinatorics of the 
sources. The combinatorial problem is to connect the dots on the upper and lower side (left- and 
right moving sources) pairwiso. 
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Fig. 16. Dominant contributions are "rainbow" diagrams, where on one side (left- or rightmoving 
sources) the same sources in the amplitude and the complex conjugate are connected to each 
other. 



and becomes correlated only when one averages over the distribution of p has the 
physical interpretation that the domainant correlations in the systems are the ones 
enhanced by large logarithms of x and are already present in the wavefunction be- 
fore the collision. Because of this the calculation of multigluon correlations is in 
fact simplified in the strong field limit |60|104) 



6.2. The ridge and the negative binomial 

We can then apply this formalism to the calculati on of the probability distribution 
of the number of gluons in the glasma^I^^II^SII^. We shall assume the "AA" power 
counting of sources that are parametrically strong in but nevertheless work to 
the lowest nontrivial order in the color sources. Formally this would correspond to 
a power counting p g^~^ with a small e > 0. In this limit, as we have discussed, 
the dominant contributions to multiparticle correlations come from diagrams that 
are disconnected for fixed sources and become connected only after averaging over 
the color charge configurations. 

Working with the MV model Gaussian probability distribution 

P'^(xt)p"(xt) 



W[p] ~ exp 



d^XT^ 



(38) 



computing the correlations in the linearixed approximation is a simple combinato- 
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rial problem. Each gluon is produced from two Lipatov vertices (see Fig. 1151 left), 
one in the amplitude and the other in the complex conjugate. The combinatorial 
factor is obtained by counting the different ways of contracting the sources pairwise 
(see Fig. [TS] right). The dominant contributions are "rainbow" diagrams, Fig. [TBI 
where on the side of one of the sources a line with momentum in the amplitude is 
connected to a line with the same momentum in the complex conjugate amplitude. 
These contributions dominate because they contain a maximally infrared divergent 
integral in the momentum circulating in the lower ("non-rainbow") side of the di- 
agram. This divergence is then regulated by the transverse correlation scale of the 
problem, Qs- When integrated over the momenta of the produced gluons one ob- 
tains the factorial moments of the multiplicity, which define the whole probability 
distribution. The result of the combinatorial exercise is that the number of con- 
tributing diagrams, each with an equal contribution, is 2'^{q ^ The probability 
distribution can be expressed in terms of two parameters, the mean multiplicity n, 
and a parameter k describing the width of the distribution. The factorial moments 
ruq = {N'i) — disc, are 

m, = {q~l)\kQ' (39) 
with parameters k and n given by 

k . (^-^ - (40) 

n = /jv — Qs'5i. (41) 

Ois 

These moments define a negative binomial distribution, which has been used as a 
phenomenological observation in high energy hadron and nuclear collisions already 
for a long time 11^. In terms of the glasma flux tube picture this result has a natural 
interpretation. The transverse area of a typical flux tube is and thus there 

are Qs'^S± — N^t independent ones. Each of these radiates particl es ind ependently 
into A^c^ — 1 color states in a Bose- Einstein distribution (see e.g. USSl)^ A sum of 
k « Nft^Nc'^ — 1) independent Bose-Einstein-distributions is precisely what defines 
a negative binomial distribution. 

An important special case of this result (and one that was first obtained in- 
dependently from the general derivation of the multiplicity distribution i are 
two gluon correlations. Similarly as in Sec. |4] one can argue thate they should be 
closely related to the two particle correlations obser ved in the final state. There are 
two major ingredients in the "glasma flux tub^^^' explanation for the ridge: 

• The long range longitudinal structure is provided by the approximately 
boost invariant glasma fields. The fields are correlated at a (transverse) 
length scale 1/Qs which enables one to estimate the strength of the corre- 
lation. 

• The narrow azimuthal structure is explained by a coUimating efrectEH 
from the strong radial fiow generated in the hydrodynamical stage. 
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Fig. 17. Left: Two-gluon correlation in a numerical CYM evaluation in the MV modeJSl. Right: 
Strength of the ridge correlation from a calculation including the initial corre lations from the 
glasma fields and a blast- wave parametrization of hydrodynamical evolutioJUSl. 



It has become conventional to parametrize the two particle correlation function in 
terms of the single particle spectra and a geometrical factor related to the number 
of flux tubes as 

C2(P,q) , ^ 1 

^_^^=«,,p„,,)_. ,42) 

wherer C2 is the two-particle correlation function 

C2(P, q) ^ ^ ^2^^ ) " ( dj/p d^pT ) ( dygd^ciT ) ■ ^^^^ 

Generally K2 is expected to be a number of o rder one. It is a constant K2 ~ 0.^ (up 
to logarithms) in the dilute limit calculatiorflS^lI^, but there is no reason for it 
to be a constant in a full calculation. The scaled correlation function has recently 
been evaluatecff^ in the MV model, with the result typically being very close to 
K2 ~ Ij but with some dependence on pt^ and qt- The result of this numerical 
calculation is shown in Fig. [17] as a function of |pt — QtI and |pt + QtI- The 
experimentally observed quantity is denoted as Ap/ y^p^ci.i^f)- Assuming a blast- 
wave parametri zatio n of the hydrodynamical evolution the value at Aip = can be 
related to K2 as 

Ap dN g2(p,q) , 1. 
(A^ 0) — . — — -— — — - ) 



= ^ (7B - (44) 
as IB 

where Kn = K2/13.5 for an SU(3) gauge theory (based on the relation between the 
multiplicity and Qs^S± discussed in Sec. HI) Here 75 is the average radial boost in 

■> After correcting for incorrect constant factors in the original calculation of LkliOJ ge for a 
discussion. 
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Fig. 18. Left: Illustration of the chiral magnetic effect. In situation 1 the spins (blue or right 
arrow) of the positively charged u quarks are aligned with the external magnetic field B, and 
those of the d quarks against it. The momenta (red or left arrow) of the right handed quarks 
are in the direction of the spin and those of the left-handed quarks are opposite to it. After an 
external color field configuration with nonzero Chern-Simons number an excess of right-handed 
quarks is present. This leads to an excess of momentum in the direction of B for the u quarks and 
against B for the d's, i.e. an electric dipole moment. Right: Correlation data showing a signal for 
a like side correlation for same charge pairs from-'^'^-l^ . The solid (AuAu) and dashed (CuCu) lines 
represent the estimates of the effect of 3-particle correlations based on HIJING. 



the framework of a blast wave model. From the RHIC data ^ one can estimate that 



Ap/v/pl^(A(p = 0) = l/y^27rcr2, with 0.64. Combining this with Eq. (|44| . 

one obtains 

for as = 0.5 and where the superscript denotes that this is a crude estimate ex- 
tracted from experiment using a blast-wave parametrization. For an average blast 
wave radial velocity Vr — 0.6, this gives ^ 1.5; for Vr — 0.7, one obtains 

^ 2. A comparison of this parametrization to experimental data is sho wn in 
Fig. [171 For more detailed phenomenological studies we refer the reader tcA^^OJ^J^. 
In spite of the increasingly detailed modeling of the final stage effects there remains 
some uncertainty in the strength of the correlation. More detailed understanding of 
the strength of the initial correlations from JIMWLK evolution and further com- 
parisons of realistic hydrodynamical calculations are still needed to clarify these 
issues. 



6.3. Chiral magnetic effect 

There have been several discussioniH^I Qf possibility of observing event-by 
event CP-violating fluctuations in heavy ion collisions. The discussion originally 
focused on the possibility of generating these fluctuations in the thermal phase 
of the collision process. However, the large longitudinal classical color f ields in 
the glasma provide a very natural framework for the phenomenon to occuJ^^^HI], 
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Parity violation in a large classical field configurations can occur when one has a 
large Chern-Simons charge Nqs: which is the time component of the topological 
current K^"; d^K^" = F^^F^^, where F^" = ^ef^'P'^Fp^. In terms of the gauge 
potentials this Chern-Simons current is 

In the strict boost invariant case Ng fi is s uppressed by topologial reasons, which 
has been seen in numerical simulation J ^ But in the presence of fluctuations that 
break boost invariance the natural (unsuppressed by the weak coupling) strength 
of Nqs is expected to be seen, although this has not yet been confirmed by an 
explicit calculation. When quarks are produced in the color field, a nonzero value 
of the Chern-Simons charge leads to an excess of quarks of a given chirality being 
produced. Due to symmetry the average of N^g over a large number of events 
will be zero, but the fluctuations will be large and measurable. The magnitude of 
the fluctuations can again be estimated using the fact that the systems consists of 
correlated areas of size l/Qs^ in the transverse plane. The number of these domains, 
each having charge ^ 1, is ^ nR^Qs^. The sum of these independent charges will 
fluctuate around zero with mag nitude ^/{N^s) ~ ~ QsRa- 

The experimental manifestation of p arity violation in heavy ion collisions hap- 
pens through the "chiral magnetic effect'^tn^ (see Fig.fTS]). In a noncentral collision 
there is an external magnetic field perpendicular to the reaction plane caused by 
the positive charges of the ions. The spins of the quarks become correlated with 
this magnetic field. In the presence of a net chirality this causes the momenta of 
the quarks to be correlated with their electric charge times the magnetic field. 
This means that one is generating, event by event, a net electric dipole moment (a 
vector) that is correlated with the reaction plane (pseudovector), a classical signal 
of CP violation analogous to the elusive neutron electric dipole moment. Such a 
correlation between the electrical charges and momenta of the produced particles 
with the reaction plane has now, after being p redic ted by the theory, been obseved 
by the STAR experimentiflQ-^EI (gee howeveJH^ for a discussion of other effects 
leading to similar experimental signatures). 

6.4. Rapidity dependence of two gluon correlation 

The same power counting argument that qualitatively distinguishes AA collisions 
from ones involving a dilute probe also applies to correlations between gluons at 
larger rapididy separations. When producing two gluons from the same BFKL lad- 
der (which is the dominant mechanism in pp and pA) momentum conservation 
causes the dominant long range (in rapidity) correlation to be back-to-back in az- 
imuthal angle. A long range correlation on the near side that leads to the ridge 
comes from diagrams that are disconnected for fixed sources and only become con- 
nected at the level of the source distribution. To calculate the rapidity dependence 
of the ridge one must therefore go beyond the leading order boost invariant clas- 
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Fig. 19. Left: Comparison of a two-particle correlation computed using Eq. I I53I1 supplemented 
with a short-range correlation contribution from PYTHIA with PHOBOS data. Right: comp arison 
Left: Rapidity correlation at LHC energies fc-p-factorization approximation. Plots frortiU^. 



sical fields. The first step in this direction is to compute the leading logarithmic 
corrections, more precisely to resum corrections that are proportional to as{yp~yq), 
where yp and yq are the rapidities of the two produced gluons. In the dense-dense 
case when the correlations are dominated by those present in the color sources, 
this is achieved by inclu ding the JIMWLK evolution between the rapidities of the 
produced gluonJU^m^. 

The crucial aspect of JIMWLK evolution for unequal rapidity correlations is 
to realize that the probability distribution W[p] should really be understood as a 
distribution for trajectories of the color charge densities (or in practice the Wilson 
lines derived from them) along path to larger rapidities. In this sense the formal- 
ism contains more physical information than just the equation for Wy[C/(xT)]; the 
probability distribution at a single rapidity y. Namely, it also gives the combined 
probability distribution for Wilson lines at different rapidities 

Wy,...y„[Uii^T),.-.,Ur,{^T)l (47) 

or in a formulation in that is continuous in rapidity W[U{x.T,y)]. We can formally 
return from the distribution of trajectories to a distribution of Wilson lines at one 
individual rapidity as 

Wypii^T)]^ j [DU{y,^T)]W[U{y,^T)]S[U{-KT)-U{y,^T)]- (48) 

Knowing the general (multiple rapidity) probability distribution will enable us to 
compute the correlations between Wilson lines, and consequently of physical observ- 
ables such as the multiplicities, at different rapidities. A transparent interpret ation 
of this can be obtained in the Langevin formulation of high energy evolutiorP^Ql^ 
where one explicitly constructs the ensemble of Wilson lines from solutions of a 
nonlinear Langevin equation. The Langevin formu lation is also the one used in the 
(numerical) full solution of the JIMWLK equatiorffSHEH 

Knowing that the correlation follows from a Langevin equation imposes an ad- 
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ditional structure (of a Markovian process) on the joint probabihty distribution: 

Wy,,y, [f/^ C/^] = G,^-,, [C/^ U'^] Wy^ [C/«] , (49) 

where the JIMWLK Green's function G is determined by the initial condition 

hm G,^_,,[C/P,[/«] = 5(t/^(xT)-C/^(xT)) (50) 

and the requirement that it must satisfy the JIMWLK equation 

dypGy^-yq [U^, t/'] = n (C/^(Xt)) Gy^^y^ [U^ , [/^] . (51) 

This JIMWLK Green's function contains ah the information, at the leading log 
level, of long range rapidity correlations in gluon production. This structure follows 
from the computation of the leading log part of 1-loop corrections to a wide class 
of observables that can be expressed in terms of correlators of the gluon fields at 
T = (or equivalently in terms of the Wilson lines) 

(0>,,„, = / [DUi{y,^T)][DU2{y,^T)]W[Ui{y,^T)]W[U2{y,^T)]0,^. (52) 

A full (numerical) calculation of rapidity correlations in JIMWLK evolu tion 
is yet to be done. The first phenomenological applications of this frameworliU^ 
follow the strategy of using the approximate linearized solution of the equations of 
motion used also in the boost invariant case reviewed in Sec. 16.21 One solves the 
two gluon correlation in terms of pp-correlatorf0, including unequal rapidity ones. 
Recognizing that due to the Markovian nature of the evolution (xt)^^ (yt)) — 
(py{xT)py{yT)), where y is the rapidity that is earlier in the evolution of nucleus 
A than y' , one can express the result in terms of equal-rapidity correlations of 
two p's. To get a realistic estimate of the resulting rapidity dependence one then 
replaces the equal rapidity pp correlators by solutions of the (running coupling) BK 
equation, denoted below as $(y,kT). The resulting fc^-factorized approximation is 

^2(p,q) = ,4 1 2 J / d'kT (53) 



IGttIO 4 pT2q2, 

iVp, ^t)^A2 iVp, Pt - kr) iVq, qr + kr) + (y?, Qt - kr 

+ *A2(y9'kT)$Ai(yp,PT - kr) $Ai(y9,qT + kT) + $Ai(y9,qT - kr) 

Although the required symmetrizations make this formula look a bit awkward, 
the structure is in fact quite simple. The double inclusive gluon multiplicity is 
proportional to four unintegrated distributions (the single inclusive one being given 
by a convolution of two). This is a demonstration of our power counting argument; 
when the distributions are assumed to be large $ ^ l/cus, this contribution to the 
correlation is ~ l/ctg. Producing the gluons from the same one BFKL ladder ("pp" 



^ We are here assuming a Gaussian distribution so that no nontrivial higher cumulants contribute. 



March 9, 2010 ijmpereview 



Small X physics and RHIC data 39 



contribution) would be parametrically a^^^ ^ 1. A result of the Markovian nature 
of the correlation is that the rapidity arguments are not symmetric; for example in 
the first term three out of the four distributions are evaluated at the rapidity yp 
(earlier in the evolution of nucleus 1) and only one at the rapidity yq. Implementing 
Eq. (|53p leads to a very weak rapidity dependence and a normalization of the 
correlation in broad agreement with the experimental results, as seen in Fig. 1191 
One could assume that the decorrelation with rapidity could be faster in the full 
nonlinear JIMWLK evolution than the linearized version assumed in deriving the 
fep-factorized approximation of the correlation. This would require combining the 
numerical JIMWLK evolution with a full GYM calculation of the multiplicity, which 
has not yet been performed. 



7. Conclusions 

We have emplasized in this short review the potential of the CGC framework as a 
description of the initial conditions of a heavy ion collision, based on weak coupling 
QCD. The CGC picture enables one to analyze DIS at small x, proton-nucleus 
collisions and the initial conditions in AA collisions in a single unified framework. 
As a weak-coupling first-principles QCD calculation there is a systematic way to 
compute higher order corrections to the results. We have seen that the basic features 
of bulk particle production in AA collisions, such as the collision energy, rapidity 
and centrality dependence can be well understood in terms of a single dominant 
momentum scale Qs present in the high energy nuclear wavefunction. An area 
of recent and ongoing activity are multiparticle correlations, which have a high 
potential of giving direct experimental access to initial state of the collision. 
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